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1 Introduction 



In this paper we are interested in the study of the following stationary 
Schrodinger system of equations: find an infinite sequence {Xm,Um)m>i and 
a potential V satisfying 

' —AUm + VoUrn + VUm = ^mUm in ^ (la) 
< °° 

-AV = Pm\um\'^ in O (lb) 

^ m=l 

where, in addition, we require that 

UmeH^in), forallm>l, F € i?o(0), (2a) 
(wm)m>i is a Hilbert basis for L'^{n). (2b) 

Here we assume that the external potential Vq and the positive numbers 
{Pm)m>i are given. 

The system (|la p - ()lbp appears in the modeling of nanoscale semiconduc- 
tor devices as part of the so-called "quantum-kinetic subband model" , which 
itself is a simplified model of the full evolution 3D-Schr6dinger-Poisson equa- 
tion N. Ben Abdallah &; F. Mehats [H [5]. In order to improve the cost of 
numerical simulation of the evolution 3D-Schr6dinger-Poisson equation and 
taking advantage of the extreme confinement of the electrons in one direction 
transverse to the transport directions, one can perform a block diagonali- 
sation of the electron Hamiltonian, thanks to a separation of the confine- 
ment and transport directions. This reduction process leads to replace the 
3D-Schrodinger-Poisson equation by a system of ID stationary Schrodinger 
equations (for the confinement) coupled to a 2D equation (for the transport). 
In that reduced model, the system of ID stationary Schrodinger equations 
is nothing but (fTa|) - (|lb|) where the parameters (pm)m>i are the sequence of 
occupation numbers, which may depend on time and space, and which are 
given by the above mentioned 2D equation (for the transport). 

On the other hand, as it is outlined in P. Zweifel ^12j, the Schrodinger- 
Poisson system of equations derives from a quantum transport equation, 
the Poisson-Wigner system. After, performing the Wigner transform to the 
former system one ends up with the following system of evolution equations 

i dttpm = -Alpm + Volpm + V^prn (3a) 

^prn{0,x) = Tpomix) (3b) 
* oo 

-AV =Y,Pm\4'm\^ (3c) 
^ m=l 
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with the condition {'4'Om\ipok) = Skm for any k,m>l. Seeking standing- 
wave solutions of the form 

^pm{t,x) := e-'^^^u^m{x) 

leads to the equations (fTa|) - (]lbp . 

Existence of solutions to the system p^ - ()lbp has been proved in R. 
Illner, O. Kavian, H. Lange [2] in the case when pm = for any m >2, and 
the same method extends to the case when pm = for any m > M for a 
given M > 2. On the other hand, a different but similar eigenvalue problem 
has been considered by F. Nier [3 [HI [7]. 

In order to solve this system of equations, taking into account the fact 
that the family {um)m>i niust be contained in ffg(Q) and, at the same time, 
has to be a Hilbert basis of L^(J7), we observe the following. Let us consider 
a fixed Hilbert basis of L^(J7), denoted by {em)m>i, such that em G Hq{Q) 
for all m > 1. For instance such a basis may be given by the eigenfunctions 
of the Laplace operator on Hq{Q), that is a family satisfying 

-Aem = IJ.mem, Cm ^ HI{Q), / ei{x) em{x) dx = 5im , (4) 

Jo. 

where the sequence of eigenvalues of the Laplace operator, with Dirichlet 
boundary conditions, is denoted by (//m)m>i- Now, saying that the family 
{um)m>i satisfies condition (l2bp means that the linear operator U acting on 
L^(ri) and defined by 

for aU m > 1, Uem--=Um, U f := '^{f\em)um. for / e L^(il), 

m>l 

(5) 

is a unitary operator, that is U*U = UU* = I. Therefore determining the 
whole family {um)m>i satisfying equations p^ - ()2bp is equivalent to find a 
linear operator U defined on L'^{Q) verifying 

U*U = UU* = I, Uej € H^{n) for j > 1, (6) 

and such that the family Uj := Uej is the family of normalized eigenfunctions 
of (fTal) where V is given by (llbp . 

In this paper we give a variational formulation of the system (jlaP "()2b p . 
yielding a solution in terms of critical points of a real valued functional 
defined on a subset of the group of unitary operators. More precisely we 
define a subset S of unitary operators on L^(0) as 

S:= {U : L^{n) ^ L^{n) ; U satisfies ([6])}, (7) 
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and then we define a functional J on S by setting 

J{U) -y^p^ [ \Vum{x)\^dx + l [ \VV[U]{x)fdx, (8) 

where Um '■= and V[U] is the solution of (fTb|) . The main purpose of 
this paper is to show that critical points of J on S yield solutions of the 
Schrodinger-Poisson system and that the minimum of J is achieved on S. 

It is clear that in order to define the potential V[U] and the functional 
J on the manifold S some conditions must be imposed on the sequence 

(,Pm)m>l- 

Our main result concerning the system of equations (jlap -()2b p is the 
following: 

Theorem 1 Let U C be a bounded domain and N < 3. Assume that 
the sequence {pm)m>i satisfies 

Pm>0, ^^^^ P"" < 

m>l 

and 

for some po > N/2 (or po = 1 if N = 1). Then the Schrodinger-Poisson 
system of equations ila \) - [2b\) has a solution obtained as the minimum of 
the functional J defined in 

In order to give a clear exposition of our global approach to the deter- 
mination of a system of eigenvectors in terms of unitary operators U, in this 
introduction we give an outline of our approach, getting rid of technicalities 
inherent to an infinite dimensional Hilbert space and to noninear problems. 
Thus, in a first step, assume that we are given a finite dimensional (complex) 

Hilbert space H of dimension n > 2, its scalar product being denoted by 
(•|-). A : H — > H is a self-adjoint, nonngeative operator (matrix), our 
aim is to define a procedure in which all the eigenvectors of A are deter- 
mined at once, to compare with a step by step construction of eigenvalues 
(and eigenvectors) through the construction of critical values of the Rayleigh 
quotient 

{Au\u) 
(u\u) 
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by a min-max proceduure. To this end, ioi ui, ... ,Un & H let us define the 
functional 

n 

F{ui,...,Un) := '^pj{Auj\uj), 

i=i 

where, as above, we assume that the coefficients pj verify pj > 0, and also, 
for the sake of simplicity of exposition (see below for the general case) , here 
assume moreover that 

Pi ^ Pj for i 7^ j. 
Define the subset (or manifold) S C by 

S := { {ui, . . . ,Un) € H"- ; {ui\uj) = 5ij } 

We claim that upon maximizing or minimizing F on the manifold S C H^, 
all the eigenvectors of A can be determined (as a matter of fact, any critical 
point of F yields such a result). 

We begin by observing that an orthonormal basis ei , . . . , e„ of if being 
given once and for all, the manifold S can be identified with the set of unitary 
matrices U such that U*U = UU* = I, where / is the identity operator on 
H: indeed it is enough to see Uj as the j-th column of U, that is to set 
Uj := Uej. Then we see that {Auj\uj) = {U*AUej\ej), and denoting by 
D the diagonal matrix D := diag(pi, . . . , pn), that is the matrix defined by 
Dcj = Pjej, we check easily that 

F{ui, ...,Un) = ^pj{Auj\uj) = ^pj{U*AUej\ej) = ti{DU*AU), 

where ti{B) denotes the trace of the operator (or matrix) B. Finally, con- 
sidering for instance the minimization of F, this can be reformulated in the 
following way: 

minimize J{U) := tT{DU* AU) under the constraint U*U = I. 

Clearly J is C°° (in fact analytic) and positive on the set 

SL(n) ■.= {U -.H ^H; U*U = I}, 

which is a smooth and compact manifold: therefore J achieves its minimum 
at some point Uq G SL(n). Now we have to show that the vectors Uj := UoBj 
are indeed the eigenvectors of A. 
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Let M : H — > H be skew-adjoint (that is M* = —M) and consider 

the one parameter group U{t) = cxp{tM)Uo for t G M; note that since 
M* = —M, one has exp(tM)* = cxp(— tM) and thus one checks easily that 
having UqUo = I, then for alH G M one has U{t) G SL(n), and consequently 
J{Uo) < J{U{t)) for all t G M. Now 

J(U(t)) = tr(DmeM-tM)AeMtM)Uo) and '^'^^^^^^^ = q, 

dt |t=o 

so that, after a straightforward calculation, we obtain 

J for all Msuch that M* = ~M, we have 
[tT{DU^MAUo) = ti{DU^AMUo). 

Setting 

B := UqDU^ , 

and using the fact that for two given matrices K, L we have tr (i^L) = 
tr{LK), we observe that tr{DU^MAUo) = tr{MAB), and that 

ti{DU^AMUo) = ti{BAM) = tr{MBA). 

Summing up, we conclude that 

for all Msuch that M* = -M, we have tr(M(^B - BA)) = 0. 

Taking M := {AB - BA)*, we conclude that BA = AB, that is 

UqDU^A = AUqDU^. 

Applying this equality to the vector Uj := U^ej, and taking into account the 
definition of the diagonal operator D, we obtain (recall that UqUq = I) 

{UQDUQ)Auj = AUoDUQUj = AUoDcj = pjAUoej = pjAuj, 

that is {UoDUQ)Auj = PjAuj, which means that Auj is an eigenvector of 
UqDUq. This implies that D{U^Auj) = pj{UQAuj), and wc sec that UqAuj 
is an eigenvector of D for the eigenvalue pj, which is a simple eigenvalue of 
D, corresponding to the eigenvector ej. This means that there exists Xj G C 
such that UqAuj = Xjcj, that is 

Auj = XjUj . 

As a matter of fact one sees that Xj G M, while Uj is an eigenvector of A 
and Uq is a diagonalization operator for A, which consists in the matrix 
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whose columns are the eigenvectors Uj. ActuaUy this procedure allows us 
to construct all the eigenvectors of A through the minimization of a unique 
functional defined on the group SL(n). Also, since F{Uo) = Ylj>i ^jPj^ o^^^ 
easily sees that different choices in ordering the numbers pj yield different 
ordering of eigenvalues and eigenvectors of A: for instance one may check 
that if the /jj's are decreasing, that is if pj > pj+i for 1 < j < n — 1, 
then one obtains the eigenvalues of ^ in a non decreasing order, that is 
Ai < A2 < • • • < An- While if pj < pj+i for 1 < j < n — 1, then one obtains 
the eigenvalues in a non increasing order, that is Ai > A2 > • • • > A„. (Had 
we began by maximizing J, the conclusion would be somehow reversed but 
analogous: if the pj''s are decreasing then the Aj's would be non increasing). 

In the next section of this paper we will show that, for a certain class 
of self-adjoint operators A, the eigenvectors and eigenvalues of A can be 
obtained through the minimization of the functional 

Jo{U) := tr{DU*AU) 

on an appropriate subset of unitary opeartors U : this is precisely stated and 
proved in section [2j In section [3] we gather a certain number of preliminary 
results used in section HI after stating the assumptions on the domain Q and 
on the sequence {pm)rm we prove theorem[Tl as well as slightly more general 
variants of the Schrodinger-Poisson systems (see theorem [10] in section [3]) . 
In section [5] we shall discuss some generalizations and state a few remarks 
about the results presented here. 

2 Global determination of eigenvectors and eigen- 
values 

In this section we consider an infinite dimensional, separable, complex Hilbert 
space whose scalar product is denoted by (•!•) and its norm by || • ||. We shall 
make the following assumptions: 

Hypothesis 1 We assume that {A,D{A)) is a densely defined, selfadjoint 
positive operator acting on the Hilbert space H , and that the domain D{A) 
equipped with its graph norm is compactly imbedded in H , so that A has a 
compact resolvent and A possesses a sequence of eigenvalues {p,j)j>i such 
that < Pj < Pj+i, each eigenvalue having finite multiplicity ruj > 1, and 
Pj — 7- +00 as j — )• 00, H being infinite dimensional. 
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Denoting by D{A^^'^) the domain of A^^"^, that is the subspace of H obtained 
upon the completion of D{A) with the scalar product (n, v) i-^ {u\v) + {Au\v), 
we recall that D{A^^'^) is dense in H. Hence we can introduce the next 
assumption: 

Hypothesis 2 We consider a fixed Hilbert basis of H , denoted by {ej)j>i, 
such that Cj € D{A^/'^) for each j > 1. 

With the Hilbert basis {ej)j>i given by hypothesis [H we consider a sequence 
{pj)j>i of real numbers such that 

Pj > 0' ^Pj l|eillD(AV2) < oo, (10) 

and we denote by D the diagonal operator defined by 

Dcj := PjCj, for j > 1. (11) 

Note that since H is infinite dimensional and A has a compact resolvent, 
while Cj G D{A^/'^), we have ||ej |[£)(-^i/2) — > oo as j ^ oo. Indeed, other- 
wise, the sequence {ej)j would be bounded in D{A^/'^), and the imbedding 
D{A^^'^) C H being compact, one would extract a subsequence {ej^)k>i such 
that ejf, f in D{A^/'^) and ej^, — > / strongly in H; in particular ||/|| = 1, 
since (ej)j is a Hilbert basis of H. But we have also ej ^ in H, and thus 
we should have / = 0. This contradiction shows that {ej)j cannot contain 
any bounded sequence in D{A^^'^). As a consequence we have pj — )• and 
D is a compact operator. 

Next we shall consider unitary operators U : H — > H which satisfy the 
following condition (this expresses the fact that the operator DU*AU is of 
trace class, see M. Reed & B. Simon [TO], volume 1, section VI. 6) 

U*U = UU* = I, Ucj G Z?(^i/2) for j > 1, 
E,>i/o,(f/MC/e,|e,)<oo. ^^^^ 

and we define the set § through 

§:={[/: F ^ F ; [/satisfies ([H}. (13) 

Remark 2 Let us point out that such operators U exist, that is § is not 
empty: indeed, for any A > 0, the operator Ux (the so-called Cayley trans- 
form of \A, see K. Yosida fU^) defined by 

Ux := {I + iXA){I - iXA)-^ 
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is a bounded operator on H and one checks easily that 

Ul = {I -\\A){I + \\AY^, 

so that UlPx = I. Moreover, for any f G D{A^/'^) we have {I + iXA)~^f G 
D(A3/2), and thus Uxf G D{A^/^). As a matter of fact, not only Ux is a 
unitary operator on H, hut one has also ||f^A/|lD{Ai/2) = II/IId(A1/2)- There- 
fore, since the sequence {pj)j satisfies [70j) . one sees that Ux satisfies [W(j 
and Ux G S. 

For a unitary operator U : H — H satisfying (fT2|) . we define Jq{U) by 

Jo(f/) := iT{DU*AU) ■.= ^pj{U*AUej\ej) (14) 

i>i 

The following result concerns eigenvectors of A: 

Theorem 3 Assume that the hypotheses\^ and\^ as well as condition 11 Op 
are satisfied. Then the functional Jq defined in |i^p achieves its minimum 
on S defined by il3\) . Then there exists Uq (zS such that 

Jo(C/o) =minJo(y), 

and Uq is a diagonalization operator for A; more precisely, for each j > I, 
the vector if j := UqCj is an eigenvector of A corresponding to the eigenvalue 
Xj := (Aifjlipj). 

We split the proof of this result into a couple of lemmas. 

Lemma 4 The functional Jq achieves its minimum on S at a certain Uq G 
S. 

Proof. Indeed consider the infimum a := inf[/g§ Jo(C/). Since S / 0, we 
have < Q < oo. Consider a minimizing sequence (C/n)„>i G S, such that 
for instance a < J{Un) < a + 1/n. Then for each fixed j > 1, setting 
u'j := Unej, we have for all n > 1 

KllV/^) = i + (^^iK)<i + ^- 

Thus, since the inclusion D{A^/'^) C H is compact, upon extracting subse- 
quences through Cantor's diagonal scheme, and denoting again this diagonal 



9 



subsequence by {Uj)n, we may assume that for all j > 1 there exist a family 
{uj)j such that for j > 1 fixed 



Uj Uj 



weakly in D{A^^'^), Uj strongly in H 



as n ^ oo. Setting UoCj := uj, one checks easily that Uq can be extended by 
linearity to the subspace spanjej ; j > 1}, and that for / G spanjej ; j > 1} 
we have 

\\Uof\\'= lim ||C/„/||2= ||/||2. 

n— >cxD 

In other words Uq is a unitary operator on (the algebraic) spanjej ; j > 1}, 
and therefore can be extended as such to the whole space H. Since for any 
m > 1 we have 

m m 

EPjiAuAuj) < liminf > pj{Au'^\u'^) < liminf Jo([/n) = a, 

upon letting m — )■ oo we conclude that Jo(f/o) < a. Thus, having UqUq = I 
and UoCj G D{A^/'^) for all j > 1, and Jo{Uq) < oo, we have C/q G § and 
MUo) = a. □ 

Next we show that Uq, given by lemma H] is a diagonalization operator 
for A. 

Lemma 5 Under the assumptions of theorem\^ let Uq be given by lemma^ 
and set Uj := Uoej, for j > 1. 

(i) Assume that k > 1 is such that 

Pi / Pk for £^ k. (15) 
Then there exist G such that Auj. = \kUk- 

(ii) Assume that k > 1 is such that for some m >2 



Pk = Pk+e for < i < m - 1, 

Pki^ Pn for n ^ {k + £ ; < i < m - 1} 



(16) 



Then there exists a unitary transformation Uk of the m- dimensional 
space Hk := span{Uoek+£ ; < ^ < m — 1} such that if 

Uk+i ■= UkUoCk+e, 

then there exists Xk+e G ^+ such that Auk+e = Xk+t'^k+t for < ^ < 
m — 1. 
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Proof. First let M : H — > H he a bounded skewadjoint operator such that 
M : D{A^/'^) — 7- D{A^/'^) is also bounded. Indeed such operators do exist 
(consider for instance i(/ + XA)~^ for A > 0). Setting U{t) := exjp{—tM)Uo 
for i € M, one checks easily that, since M* = —M, one has U{t) € § for all 
t, and thus the function g{t) := Jo(C/(t)) is well defined, is of class and 
achieves its minimum at t = 0. However since 

g{t) = ti:{DU^ exp{tM)Aexp{-tM)Uo), 

one concludes that 

g'{0) = tr{DUoMAUo) - tT{DUoAMUo) = (17) 

for all bounded operators M : H — > H such that M* = —M and M is 
also bounded from D{A^/'^) into itself. In the same way, if we consider a 
bounded operator L : H — > H such that L = L* and L is also bounded 
from D(A^^'^) into itself, upon setting M := iL, we conclude that ()17p yields 

tr{DUoLAUo) = tT{DUoALUo), (18) 

for all such operators L. 

Note that the above relation (flT]) yields that 

Y,PjiUoMAUoej\ej) - /5,(C/*^MC/oe,|e,) = 
i>i j>i 
Y^Pj{MAuj\uj) -Y,Pji^Muj\uj) = 0, 
i>i i>i 

that is, since M* = —M, 

-^Pj{Auj\Muj) - ^pj{AMuj\uj) = Re^pj{Auj\Muj) = 0. 
i>i j>i j>i 

(19) 

Analogously using (jlSp one obtains in the same way 

""^^ Pj{Auj\Luj) = ""^^ Pj{ALuj\uj) <^=^ Im pj (Auj \Luj) = . (20) 
j>i i>i j>i 

At this point, in a first step, assume that the integer k is such that 
condition (fTSj) is fulfilled. Consider an integer n ^ k, so that Pn ^ Pk, and 
define the operators M and L in the following way 

Muk := n„, Mn„ := -Uk, 
< Lufc := Un, Lun := Ufc, (21) 
^Luj = Muj = for j {A;, n}. 
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Clearly M and L satisfy the required conditions above, and using (|19p. 
with our choice of the operator M, we get (p„ — pfc)Re(Anfc|un) = 0, that 
is, since Pfc - Pn / 0, 

Re(^tifc|n„) = 0. 

Upon using ()20p . with our above choice of the operator L and the fact that 
Pn — Pk 0, analogously we have that 

lui{Auk\un) = 0. 

So, from the above two relations, we infer that {Auk\un) = for all n such 
that Pn Pk-, that is 

Auk € span{ti„ ; n ^ k}^ = spanjufc}, 

where we use the fact that the family {uj)j is a Hilbert basis of H, being the 
image of the Hilbert basis {ej)j under the unitary operator Uq. This means 
that Auk = ^kUk for some G C, but since A is a nonnegative self-adjoint 
operator, as a matter of fact we have Afc > 0. 

Next assume that the integer k is such that the coefficient pk has mul- 
tiplicity m > 2, that is condition (I16p is satisfied. Arguing as above, we 
consider the following operators M and L: for n ^ {k + j ; < j < m — 1} 
and < £ < (m — 1) fixed, set 

Muk+e ■■=Un, Mun := -Uk+e, 
< Luk+e ■■= Un, Lun := Uk+e (22) 
^Luj = Muj = for all j ^ {n,k + i}. 

Then, proceeding as above, we conclude that {Auk+e\un) = for all n 
{k + j ] < j < m — 1}, that is: 

Auk+e G (span {un ] n ^ k + j, <j <m - 1})^ 

that is 

Auk+e € span {uk+i ; < i < m — 1} . 

This means that if we set Hk := span {uk+i ; < i < m — 1}, then A : 
Hk — > Hk is a self-adjoint operator on the finite dimensional space Hk- 
Therefore there exists a unitary operator Uk, acting on this space, such 
that if for < £ < m — 1 we set Uk+i = UkUk+e = UkUoCk+i, we have 
Auk+e = Xk+eUk+e for some Xk+e > 0. □ 
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As we may see from the above analysis, when all the pj's are distinct, 
then Uq, any unitary operator which minimizes Jq, is a diagonalization op- 
erator for A. However in the general case, when some of the coefficients pk 
have multiplicity > 2, it is possible that one has to impose a unitary 
transformation in the space 

Hk := span{C/oeA,.+<? ; < i < nik - 1} 

in order to have the operator A diagonalized. In other words, one may find 
a unitary operator Uk on Hk such that if Ak := A^h^. is the trace of A on H^, 
the operator U^A^Uk is diagonal. Thus since pk = Pk+i for < ^ < m — 1, 
if we denote by U the unitary operator obtained through the composition 
of all such operators Uk and C/q, one has Jo{U) = Jq{Uo). More precisely, 
we can we state the following corollary, which ends the proof of theorem [S) 

Corollary 6 Under the assumptions oftheorem\^ let Uk he given by lemma\M 
when k > 1 is such that il6\} is satisfied. Then the operator Uq defined by 
Uoek = Uoek when k satisfies [T5\} . and 

UoCk+e := UkUoCk+e, for < £ < m — 1, when < f7g|) is satisfied, 

belongs to S, while Jo(f^o) = Jo{Uo) 'md UqAUq is diagonal. Setting ipj := 
UoCj for j > 1, then there exists Xj > such that Aipj = Xj^pj. 



3 Preliminary results for Schrodinger— Poisson sys- 
tem 

In this section we prove an existence result regarding the system of equations 
([Ia])-(l2b]). We shall assume that 

Q C is a bounded domain and that < 3, (23) 

and we endow the (complex) space L^(r2) with its scalar product denoted 
by (-j-) and its norm || • |[. Let Vq be a real valued potential such that 

^0+ E L^{n), Vq- £ LP°{n) for some Po > y and po > 1- (24) 

Then we define an unbounded operator (A, D{A)) by setting 

D{A) := |n € H^{n) ; -An + Vqu G L^{n)^ , Au := -Au + Vqu. 

(25) 
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This operator is self-adjoint, has a compact resolvent, and there exists a 
Hilbert basis of eigensystem denoted by {Xm,^m)m>i, that is (here 6mn 
being the Kronecker symbol) 

-Aipm + Vo'Pm = >^m'Pm, (fim & H^in), / ipm{x)(pn{x) dx = 6mn- 

It is well-known that by Weyl's theorem there exist two positive constants 
ci, C2, depending on Q, and Vq, such that for all integers m > 1 one has 

(See for instance [3], chapter 5, § 3, where the case of Neumann boundary 
conditions is also treated). For this reason, as far as the sequence {pm)m>i 
is concerned, in order to ensure the finiteness of the functionals we are going 
to minimize, we assume that the growth condition Q is satisfied. 

For a given unitary operator U : L'^{Q) — > L'^{Q,) we shall denote by 
V := V[U] the potential defined by the Poisson equation (here \Uipj \ denotes 
the modulus of the function Uipj) 

- AV = ^Pj\Uipj\'^ in n, VeH^in). (26) 

It is noteworthy to recall that, thanks to the Sobolev imbedding theorem we 
have H^{n) C L^' (Q), where 2* = 27V/ (iV - 2) when TV > 3, while 2* can be 
any finite exponent if = 2, and 2* = oo if = 1: thus, when < 3, the 
right hand side of the above equation belongs to some L'^[Q) with q > A^/2, 
and more precisely a classical regularity result (see for instance [1]) states 
that there exists a constant c > such that if y G //q (il) satisfies —AV = f 
and / G L^in), then 

II^IIOO <C|1/Ilg. (27) 

Therefore, if V := V[U] is given by ([26]) we have V € L°°(0) when N <3 
(see below lemma [9]) . Note also that by the maximum principle we have 
V[U] >Omn. 

It is useful to consider the Sobolev space Hi endowed with the norm 
II ■ IIhi: 

Ml := |u G H^{n) ■ \\u\\'^^ := \\Vuf + j ix)\u{x)\'^dx < oo 

The imbedding Hi C L^(0) is compact. Note that since the eigenfunctions 
ipm belong to L°°(il), we have ipm G Hi. 
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Regarding the manifold S defined in d?]), we have to modify it shghtly, 
as we did in section § 2. More precisely we shall consider unitary operators 
U : L2(0) — > L'^{n) such that 

U*U = UU* = I, Uipj £ Ml ioi j > I, pj{U* AU if j\ip j) <oo, (28) 

and we consider the manifold defined by 

§:={[/: L'^{Q) L^{Q) ■ [/satisfies ([28])}. (29) 

We denote by D the diagonal operator acting on L^(r2) defined by D^pj = 
Pjfj, and for ?7 € § we define the functionals Jq and Ji as follows: 

Jq{U) ■.= tiiDU*AU) = S2pj f (\^U^j\'^{x) + Vo{x)\Uipj\'^{x))dx (30) 
and 

MU) -.= 1 f \VV[U]Wx)dx 
^ Jo. 

= {-\aV[UIV[U]) (31) 

where, in the last equality of (|3ip . by an abuse of notation, we denote by 
V\U] the (linear) multiplication operator / i->- y[C/]/. Since we assume 

< 3, we know that V\U\ € L°°{yt) as seen above, so that this operator 
acts boundedly on L^(0). 

Note that here the potential Vq may have a negative part, so at some 
point we will need to ensure that the functional Jq is bounded below, and 
that is coercive in some sense. More precisely we have: 

Lemma 7 There exists C > such that for any U €S one has 




Proof. Assume that = 3 (the case N < 2 can be handled in a similar 



15 



way). For t > and u G Hq{'D,) such that ||n|| = 1 we have 

/ Vi^\u\'^dx= Vf^\ufdx+ Vf^\u\'^dx 

Jn •^[Kr>t] J[Vo~<t] 

- >t] o~ \'^\'^ dx + t / l^l^da; 

< Ci(iV)meas([yo-" > t])^ \\V^\\lpo \\Vuf + t 

where we have used Holder's inequahty twice (once with N/{N — 2) and 
{N/{N - 2))' = N/2, once with po and N/2, where 9 = {2/N) - (l/po) = 
(2/3) — (1/po) > 0). Wc used also Sobolev's inequality ||n||2* < C[|Vn||. 
Now, since V^" G L*'"($7), we know that meas([FQ~ > t]) — ^ as t — ^ +oo. 
We choose t > large enough to ensure that 

Ci{N)meas{[%- > t]f HFo'lU™ < ^- 
Then we have for all u eMi 

f \Vufdx+ f Vo\u\'^dx>h\Vuf + [ V^\ufdx-t. 

Applying this to u := U(pj, multiplying by pj > and calculating the sum 
over j yields the inequality claimed by our lemma, with C := ^ X]j>i pj- D 

It is well known that the fact that the functional u ^ !|Vu|p is weakly 
sequentially lower semi-continuous (l.s.c.) on Hq(^[1) plays a crucial role in 
many, if not all, minimization problems. Regarding the functional Jq we 
need an analogous property which is stated below: 

Lemma 8 The functionalJo is "weakly sequentially lower semi-continuous" 
in the following sense: let {Un)n>i be a sequence in S such that for some 

R > and all n > 1 one has Jo{Un) < R- Then there exists a subsequence 
{Un^)k such that for any fixed j > 1 one has Un^ipj uj in Mi as k ^ +oo, 
and if we set uj'' := Unf.<^j and we define a linear operator U by setting 
Uifj := Uj we have u^^ — >■ Uj strongly in L^(J7) and 

lim /" = / V^\uj\^dx, U eS, Jo(?7) < liminf Jq (?/„,). 
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Proof. Assume = 3. Thanks to lemma [71 we know that 

/ (^[VC/„99jf + 2V(^\Un^j\^'^ dx<2R + 2C =: Ci. 

This imphes that for each j > 1 fixed the sequence (tt")n := {Un^j)n is 
bounded in Hi, more precisely 11^" < C / pj. By using Cantor's diagonal 
scheme and the compactness of the imbedding Hi C L?'{VL), we may extract 
a subsequence denoted by (u"*)fc>i such that 

^"fc ^ Uj weakly in Hi, 
ti"*" Uj strongly in 
y^Uj Uj a.e. m iZ, 

as A; — )• oo. For any m > 1 fixed, we have 

mm oo 

^PjWujWm^ < liminf^pjlln"*^!!!,^ < liminf ^pjHu"*!!^^ < C. 
j=i j=i j=i 

and finally 

oo oo 

^PjWujWm^ < liminfX^^ill^rilBii < C'^ (32) 

^ — ' fe— >oo 

Setting Uipj := Uj, one checks easily that U can be extended by linearity to 
the subspace span{ipj ; j > 1}, and that for / € spanjc^j ; j > 1} we have 

\\Uff= hm ||;7,/||2 = ||/||2. 

ra— )-oo 

In other words U is a unitary operator on (the algebraic) span{ipj ; j > 1}, 
and therefore can be extended as such to the whole space L2($7). Then ()32p 
shows that C/ G S. 

We note also that in particular we have 



7>1 



Uj{x)'^dx < oo. (33) 



Since we assume A = 3, the strong convergence of u^*" Uj in L (fl) 
implies (through Holder's inequality, or interpolation between L'^{Cl) and 
L^{Q,)) that for any fixed p < 3, and any j > 1 we have that u^'' Uj 
strongly in L'^p{Q) and a.e. in Q, and thus — > \uj\'^ strongly in L^(0). 
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Since Vq € LP°{Q) and po > 3/2, taking p ■= p'q = Po/{po — 1) we conclude 
first that 

lim / [^(ix = / Vo'lujl'^dx, 

and then thanks to (|33p and the monotone convergence theorem, 

hm Y,Pj / ^o"l^"1^^2; = / VQ-\u/dx. 

fe— >-oo ^ — ' In — In 

From this and (1321) it is clear that 



CXD „ 

V/5j / (iV-Ujf + VQ+{x)\uj\) dx < liminf Jo([/„J 
^ — ' In ^ ' fe— i>oo 

+ lim / Fo^ln^" 

7^ J<^ 



^dx 



< liminf Jo(;7„J 

fc— )-oo 



3'- 

which means that Jq{U) < liminf^^oo Jo(C/„^,), as claimed. □ 
Regarding the functional Ji we have the following result: 

Lemma 9 The functional Ji is "weakly sequentially continuous" in the fol- 
lowing sense: let {Un)n>i be a sequence in S such that for some R > and 
all n > 1 one has 

Pj\\Un^j\\M, < R, 

and such that for any fixed j > 1 one has Un^j uj in Hi as n +oo. If 
we set Uj := Un^j and we define a linear operator U by setting Uipj := Uj 
we have V[U''] V[U] strongly in L°^{Q)nH^{n) and Ji(C/") ^ JiiU) as 
n — > oo. 

Proof. Let us assume = 3, as the cases < 2 is easily handled anal- 
ogously. Since ^ Uj in Hi, we have n" Uj strongly in L^(0). By 
Sobolev imbedding theorem we conclude that for 3 < p < 6 for some con- 
stants c > and 9 := 3{p — 2)/{2p) we have 
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From this we infer first that — )• |%|^ in as n ^ oo, and next 

that for any m > 1 fixed, and setting go 1/^ > 1 and q'q = qo/{qo — l) < oo 
we have 



j>m 



<Y^p,\\u-\\l<cY^p,\\Vu-\ 
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j>m J \i>m J \j>m 

From this it follows that 

Pj\u'^\^ Pj\uj\^ strongly in LP/^{a), 

i>i j>i 

for any p € (3,6). Therefore thanks to ([27]) we have that V[Un] V[U] 
strongly in L°°(il) and finally 

[ ypn] \u]\'^dx ^YpJ [ ^[^] {ujl'^dx, 
as n ^ oo, which means that Ji{Un) — > Ji{U). □ 

4 Existence of solutions for the Schrodinger— Poisson 
system 

In this section we solve the following Schrodinger-Poisson problem: 

-Atim + VoUm + Vum = XmUm in O (34a) 

oo 

-AV = /Omi'UmP in ^ (34b) 



771=1 



and moreover 



Um G i?o(f^)' for all m > 1, V e H^i^), (35a) 
("Urn, )r?i>i is a Hilbert basis for l2(J7). (35b) 



The main result of this section is: 
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Theorem 10 Assume that the hypotheses i23\ )-(24^, as well as condition 
^ are satisfied. The functionals Jq and Ji being defined in l[30\) - [31]) . we 
set J{U) := Jq{U) + Ji{U) for U G § given by flPj) . Then J achieves its 
minimum on § and there exists C/q G S such that Jq{Uo) = mmjj^s Jq{U), 
and the family uj := Uoipj is solution to ^^3J^-(35W- Moreover if V G 
Hq{Q.) satisfies by —AV = Ylm>i Pm\um\'^, then the eigenvalues Xj are given 
by 

Xj := {-Auj + VoUj + Vuj\uj) = / \Vuj\'^dx + / (Vo + V) Uj{xfdx. 

We split the proof of this theorem into several lemmas. First we show 
that J achieves indeed its minimum. 

Lemma 11 The functional J achieves its minimum on § at a certain Uq G 
S. 

Proof. Since § ^ and Jq is bounded below (see lemma [7]), so is J and the 
infimum 

a := inf J(U) 

is finite. Consider a minimizing sequence {Un)n>i G such that for instance 
a < J{Un) < a + l/n. In particular Jo{Un) < 1 + and thanks to lemmaO 
there exists a subsequence (which denote again by Un) such that if we set 
setting n" := Un^j for each fixed j > 1, we have, for all 2 < p < 2*, 

Uj Uj weakly in Hi, 
— )• Uj strongly in LP(0), 
yu^ ^ Uj a.e. in Q, 

and the operator Uq being defined by Uoipj := uj, we have Uq £ 8> and 

Jo(C/o) <liminf Jo(f/n). 

On the other hand, thanks to lemma [H we know that ^[[/n] ^[Uo] in 
L°°{n) and that Ji([/„) ^ Ji{Uo). Thus 

a = liminf J(C/„) = liminf Jo(C/„)+ lim Ji(C/„) > MUo) + Ji{Uo) = J{Uq), 
that is, since Uq G S, we have Jo{Uo) = a = inf{/g§ J{U). □ 
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Lemma 12 Let Uq be given by lemma and let M : H — > H be a 
bounded skewadjoint operator such that M : Hi — > Hi is also bounded. Set 
U{t) := exp(— tM)C/o for t € M, and go{t) := Jo(C/(t)). Then go is of class 
and 

g'oiO) = -2ReY,PjiAUoiPj\MUo^j). (36) 
j>i 

Proof. First one checks easily that, since M* = —M, one has U{t) G § for 
ah t, and thus the function 

go{t) := MUit)) = Y,PMU{t)ipj\U{t)ipj) 

is weh defined and is of class C^. Since U'{t) := dU{t)/dt = -M exp{-tM)Uo, 
one sees that 

9'oit) = Y.Pj{AU'{t)ip,\U{t)ip,) + Y.PMU{t)ipj\U'{t)ip,), 
and finally, 

5o(0) = -^PjiAMUoipjpoipj) -J2pji^UoiPj\MUoipj), 
j>i i>i 

which yields our claim since A* = A. □ 

We have an analogous result concerning the functional Ji: before show- 
ing this, we need to show that the mapping U V[U] is smooth. 

Lemma 13 Let Uq be given by lemma and let M : H — > H be a 
bounded skewadjoint operator such that M : Mi — > Hi is also bounded. Set 
U{t) := exp{-tM)Uo fort G R. Denoting by V{t) the mapping t ^ V[U{t)], 
then t I—)- V(t) is of class from M into L°°(il) n Hq{Q) and denoting by 
W the solution of 

-AW = -2Re ^ pj {MUoVj)U^, W € H^{n), 
we have V'{0) = W. 

Proof. The fact that for all T > and t G [-T, T] we have 

||Mexp(-tM)[/o93j||Hi < ||Af ||£(]H[i) exp(r||Af||£(e^) ||[/oV9j ||hi , 
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shows that for any p G (3, 6) (see the proof of lemma [9]) the mapping 

is of class from {-T,T) — > LP/^{n) C H-^{n), and thus using (HH), 
the mapping t ^ V{t) := V[U{t)] is of class from {-T,T) into L°°{n) n 
Hq{Q). The calculation of V'{0) is straightforward. 

□ 

Now we can state the following result, which will allow us to characterize 
Uq given by lemma [TTl 

Lemma 14 Let Uq be given by lemma and let M : H — > H be a 
bounded skewadjoint operator such that M : Mi — > Hi is also bounded. Set 
U{t) := exp(-tM)[/o for t G R, and gi{t) := Ji{U{t)). Then gi is of class 
and 

g[iO) = -2ReY,P3iV[Uo]UoiPj\MUoipj). (37) 
i>i 

Proof. Thanks to lemma [13] one checks easily that the function 

51 (t) := Ji([/(t)) = ^ \VV[U{t)]fdx = ^{-AV[Uit)]\V[Um 
is C\ and that denoting by W := V'{0), we have (with Vq := V[Uo]) 
g[iO) = {-AW\Vo) = -2Re Vpj / {MUoipj)Wp'jVodx 

= -2Re^pj{MUQ^j\VQUo^j), 

i>i 

where we use the fact that Vq := V\Uq\ is real valued. □ 

The following result is analogous to lemma [SJ the only difference is that 
due to the presence of the nonlinear term we have to check that when some 
Pfc has multiplicity m > 2, we can still proceed as before. 

Lemma 15 Under the assumptions oftheorem\10i, let Uq be given by lemma [Tl\ 
and set uj := UQipj, for j > 1. Then the conclusions of lemma\5\ hold. 

Proof. The proof is very much the same as in lemma [5l so we give only the 
outline and the changes to be made. With the notations of lemmas [12] and 
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[HI we set g{t) := J{U{t)) = go{t) + gi{t). Since g{0) < g{t) for all t G M, 
we have ^'(O) = 0, that is 

ReY^Pji^UoVjlMUoifj) + Re ^ pjiV[Uo]Uoipj\MUo^j) = (38) 
j>i j>i 

for all bounded skew-adjoint operators M such that M : Hi — > Hi is also 
bounded. 

In the same way, if we consider a bounded adjoint operator L such that 
L : Hi — > Hi is also bounded, we may set M := iL and conclude that (p8]) 
yields 

Im Pji^Uo^j \LUoVj) + Im ^ pj{V[Uo]Uo^j \LUoipj) = (39) 
i>i i>i 

for all such operators L. 

At this point, in a first step, assume that the integer k is such that 
condition (fT5]) is fulfilled. Choosing M and L as in (j2T]) . and proceding 
exactly as in the proof of lemma \5\ using the fact that Pk — Pn 7^ 0, we 
conclude that 

{Auk + V[Uo\Uk\Un) = 

that is 

Aufc + y[C/o]ufc E span{n„ ; n 7^ fc}^ = spanjufc}. 
This means that 

Auk + V[Uo\uk = XkUk 

for some G C, but since A + V[Uo] is a self-adjoint operator, as a matter 
of fact we have € M. 

Next assume that the integer k is such that the coefficient pk has mul- 
tiplicity m > 2, that is condition ()16p is satisfied. Arguing as above, we 
choose the operators M and L as in (p2]) . and conclude that 

(^'Ufc+^ + l^[f/o]^ifc+^|^tn) = 

for all n {/c + J ; < j < m — 1}, that is: 

Auk+e + y [C/o]nfc+£ G (span {n„ ; n^k + j, <j <m - 1})^ . 
This means that if we set Hk ■= span{iifc-|-j ; < i < m — 1}, and 

Aqu := -Au + (Vo + y[[/o]) u 
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then Aq : Hk — > is a self-adjoint operator on the finite dimensional space 
H^. Therefore there exists a unitary operator U^, acting on this space, such 
that if for < ^ < m — 1 we set Uk+t = UkUk+i = UkUoipk+e, we have 
AoUk+i = Xk+eUk+e for some Xk+i G M. 

However, since Uk is a unitary operator on we have 

m—1 m—1 

\Uo(fk+e\'^ = \UkUo(pk+e\'^ 



/ , I ^ u ri;-r-<- 1 / J 

£=0 ^=0 



and thus 



fc+'m— 1 jn— 1 m—1 

i=k l=Q ^=0 

This means that if we set U^^j := Uqcpj if j ^ {k + i ; < i < m — 1 and 
Uj = UkUo^Pj if k < j < k + m — 1, we have F[f/o] = V[Uk], and finally this 
implies that 

□ 

As we may see from the above analysis, when all the pj's are distinct, 
then C/q, any unitary operator which minimizes Jq, yields a solution to the 
Schrodinger-Poisson system. However in the general case, when some of the 
coefficients pk have multiplicity > 2, it is possible that one has to impose 
a unitary transformation Uk in the space 

Hk := span{C/o(/?A,.+<? ; < £ < - 1} 

in order to obtain a solution (note that these unitary transformations do 
not change the value of F[C/o]). 

In other words, one may find a unitary operator Uk on Hk such that if 
Ak := A^jj^ is the trace of A on Hk, the operator U^AkUk is diagonal. Thus 
since pk = Pk+i for < £ < m — 1, if we denote by U the unitary operator 
obtained through the composition of all such operators Uk and f/o, one has 
Jq{U) = Jq{Uq). More precisely, we can we state the following corollary, 
which ends the proof of theorem [151 
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Corollary 16 Under the assumptions of theorem \1(A let Uk be given by 
lemma [73 when k > 1 is such that U6\) is satisfied. Define the operator Uq 
by Uoipk = Uocpk when k satisfies [T5\) . and 

Uoifk+t := UkUo(pk+£, for < £ < m — 1, when is satisfied. 

Then Uq belongs to S, while Jo{Uo) = Jo{Uo) and V[Uo] = V[Uo]. Moreover 
setting Uj := Uoipj for j > 1, there exists Xj G M such that 

—Auj + (Vo + V^[C/o]^ Uj = ^jUj, Uj G Hi, {uj\uk) = 5jk 

and moreover V := F[C/o] satisfies 

5 Further remarks 

The one dimensional case d = 1 is particularly simple to handle, using a 
completely different method. Indeed we point out that for a given potential 
V € C([0, 1]) the spectral sequence {Xk,^k)k>i 

-ip'l + {V + V)^k = Xk^k, (^fc(O) = (^fc(l) = 0, (^Uo)>o, 

is well defined, and ipk^fjdx = 6kj, each eigenvalue being simple. Using 
the simplicity of the eigenvalues, it is known that the mapping V ipk 
continuous fron C([0, 1]) into L^(0, 1) (see J. Poschel & E. Trubowitz [9]). 
If the coefficients {pj)j>i satisfy 

Pj > 0, y^pj =: M < oo 
i>i 

one can easily see that the mapping F : C([0, 1]) — > L^(0, 1) 

V^F{V) ■.= J2pk\Vk\^ 

k>l 

is continuous. 

Now consider the mapping B : L^(0, 1) — > L^{0, 1) defined by Bf := v 
where v G C^([0, 1]) is given by 

-v" = f in (0,1), t>(0) = t;(l) = 0. 
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Clearly B can also be considered as a linear mapping on C([0, 1]), and B is 
compact. Also observe that the second equation in (jlbp is equivalent to find 
V G C([0, 1]) such that 

V - BF{V) = 0. 

Denoting by Tiy) := BF{V), we know that / — T is a compact perturbation 
of the identity on C([0, 1]) and one can check easily that there exists R > 
such that for any E [0, 1] 

V - 9T{V) = ^ ll^lloo < R. 

Thus the invariance by homotpy of the Leray-Schauder topological degree 
implies that for all G [0, 1] we have 

deg(/ - er, B{0, R),0) = deg(/, B{0, R),0) = 1, 

which means in particular that deg(/ — T, B{0, R),0) = 1. Therefore there 
exists at least one V G C([0, 1]) such that V — T{V) = 0, that is the system 
P^ - ()lbp has at least one solution, when Q, = (0, 1). 

We point out also that in the case in which pm is a function of Am, 
for instance pm '■= exp(— A^) (see F. Nier [8]), the same approach can be 
applied. 
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